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On Sound Channel in Stratosphere 


By Yoshihito TAKESADA 


Kyoto Gakugei University 
(Read May 4, 1961; Received May 25, 1961) 


Abstract 


The sound channel will be formed in accordance with the wind current and 
temperature distribution. When the sound ray has the same direction of the 
propagation with the jet stream, the velocity minimum in the vertical distribution 
of the sound velocity exists on about 23km in height. And if the sound ray 
radiating from the ground, it is refracted at about 13km in height and then comes 
back to the ground. When the sound ray has opposite direction to wind stream, 
the velocity minimum of the vertical sound velocity distribution will occur at the 
altitude of 12 km. 


1. Introduction 


In the past, the study of sound propagation in the atmosphere has been carried 
out by S. Fujiwara, S. Shono, F.J.W. Whipple, B. Guttenberg, O. Meisser, N.A. Haskell 
and others. On the other hand, a short period phase which is travelling through the 
ocean with the velocity of sound is already known as T phase. As to the stratosphere 
“and ionosphere, from the standpoint of Acoustics, there are mutual relations between 
the meteorological phenomena and very low frequency sound. According to the data 
of Meteorological agency of Japan during 1957 and 1958, firstly, we investigated the 
most monthly probable values of the wind velocity without considering its directions 
over the upper atmosphere in Japan because of the strong wind directions are recognized 
between the angles of about 260° and 290°. Next, we investigate the general temperature 
distributions for the latitude over the Japanese islands, and then considered the effect 
of the wind velocity on the vertical distribution of the infra-sound velocity which 
depends upon the temperature only. From the gradient of the vertical distribution 
of the sound velocity, the paths of the rays will be defined. If the sound waves are 
generated in upper atmosphere the waves propagate into the channel, and there will 


be acoustic cloud. 
2. Sound Velocity and Intensity above the Stratosphere 


h 
The equation of the atmospheric pressure is given as P=P,—\ ogdh. Then 


P=P,—[pgh], where p and g are mean values in any altitude #. General formula of 


the sound velocity is given as C= a uC Introducing the density o(=75,) and above 


: : ue (ep t\ roegh t 
pressure formula to this Laplace formula we derive C= ‘a ce (1+5%3)- a : (14553), 
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where P, and p) are the pressure and the density at 0°C, 1 Atm. respectively. The distri- 
bution of the sound velocity is given by above formula. From this formula the velocity 
minimum exists at about 10km and 70km~90 km in height. If we introduce the wind 
current to this distribution of the sound velocity, the practical velocity is determined. 

An edge tone is generated in upper atmosphere by achain of mountains. For the 
purpose of propagation of the sound, the variation of the density 4o must be smaller 
tnen the original density o). The value log;)¢.=—4.3---+:- (1) at an altitude of about 
20km. Radiation energy of the sound is given as E=p ane’ tenes (2), where P is 


the pressure amplitude of the sound, x is a distance from the sound source to the 
observation point and pc is the acoustic impedance of the medium. On the other hand, 


total energy of the sound is E=+pck tenes (3)3 where é is the particle velocity of the 
medium. And as the relation between particles velocity and sound velocity, we have 
&=<CG settee (4). In this case, gradient of wave form G is equal to 4p/p,. From the 


relations (3) and (4) the total energy of the sound can be represented as E<+CG* at () 
Then from the relations (1), (2) and (5) we can derive the relation as shown in Fig. 1. 
It may be considered that the density waves can be superposed in the sound channel. 
When the value of 4o reaches about 107, the wave form of the sound will be trans- 
formed into the wind current. 


Lower Limit of Density Variation, ap 


/ 1 a 
Sound Pressure and Fluid Pressure dyne/er* 


(07 «10% ~ 7403)" o*” “yo? 


———; sound pressure. ---; fluid pressure. 
Fig. 1. Relation of the limiting values of sound pressure 
and maximum variation of the density on the axis of the 
sound ray at the distance of 100cm from the source. 


3. Sound Channel 


We can draw the figure of the long propagation distance of the infra-sound depends 
upon the distribution of the sound velocity. Sound ray is provided by Snell’s law 


in correspondence to the sound velocity distribution with altitude. If we take a given 
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division of horizontal layer in any free space of atmosphere, the following relations 


are obtained, i.e. COSHg 2 t , Yé=tan Git Oi-1 ,as shown in Fig. 2, where the notations 
COSs6;_; Cots) Se te 


of 0;, ci, x; and y; are angles of the sound ray to the horizontal line, sound velocity, 
horizontal distance and vertical height in layer of the 7 th order, respectively. The 
notation of Cy is the sound velocity on the position where the progressing sound ray 
is horizontal and sound velocity minimum is denoted as C. In this case, Snell’s law 


is given as cosé= & : 
0 


Altitude km 


Propagation Distance &m 


Fig. 2. Model of the sound channel. Direction of the sound 

propagation is opposite to the wind direction. If the critical angle 

6=33°, the horizontal propagation distance from the source to the 
ground surface is about 48.5 km. 


4. Sound Ranging in the Stratosphere 


Now, we shall explain a long distance ranging of the sound in the stratosphere. 

Over the Japan islands in winter at altitude about 13 km, the westerly current of the 
_-wind is strong as already explained by Chicago University Groups. When these wind 
velocities are added to sound velocities at each altitude respectively, we obtain the 
most possible distributions of the sound velocity as shown in Fig. 3. “According 
to the tendency of the distributions of these sound velocities, we divide into three 
zones the space up to 55km in height. As to each of these zones, we consider the 
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propagation of the sound in accordance with Snell’s law. In this sound ctleainaba 
case, if the sound source exists on the ground, a greater part of the sound will be 
refracted at the boundary between the troposphere and stratosphere. On the other 
hand, when we apply the distribution of the wind velocity in upper stratosphere to 
the distribution of the sound velocity and if the sound source exists in the stratosphere, 
the sound will propagate a long distance along the axis of the sound channel in the 
stratosphere. Accordingly, when the sound propagates in the wind direction, the paths 


of the infra-sound are given as Fig. 4. 
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Fig. 4. Ray paths of the sound propagation. Direction of the sound propagation is the same 

to the wind direction. Distribution of the sound velocity is obtained from Fig. 3. In the 

troposphere ; max. sound velocity Cy>=375 m/sec, min. sound velocity C,=340 m/sec and the critical 

angle is 25°. In the lower stratosphere ; max. sound velocity is 375m/sec, min. sound velocity is 

300 m/sec and the critical angle is 37°. In the upper stratosphere; the sound velocity depends on 

temperature and assumed wind velocity. Max. sound velocity C;=395 m/sec, min. sound velocity 
C,=300 m/sec and the critical angle is 40° 40’. 


5. Conclusion 


The altitude of the sound channel axis depends on the direction of the sound ray 
and the circulation current of the medium over the earth surface. The possible axis 
of the sound channel in stratosphere is taken to be the space between the altitudes 
of about 12km and 23km. And in the domain of the sound channel, the existence 
of comparatively weak infra-sound can be predicted. From the consideration of above- 
mentioned channels, we can estimate that there is an axis of the sound channel on about 
70km~90 km in height through a calculation of the sound velocity and the wind 
distribution over the stratosphere. 

Finally, the author wishes to express his gratitude for the discussion of Dr. Morton G. 
Wurtele in University of California. And he also wishes to express his sincere thank 
to Honorary Professor Dr. M. Hasegawa for his constant guidance. 
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Time Variation of Cosmic Ray Intensity in the Antarctic Region 
I. Atmospheric Temperature Effect 


By Taiichi KITAMURA 


Geophysical Institute, Kyoto University 
and 


Masahiro KODAMA 


The Institute of Physical and Chemical Research, Tokyo 
(Read October 30. 1960; Received June 10, 1961) 


Abstract 


Cosmic ray observations at Syowa Base, Antarctic, have been carried out by 
using a Neher type ionization chamber during 1957 and a cubical meson telescope 
since 1959, respectively. 

The atmospheric temperature effect on cosmic ray intensity thereby obtained 
has been investigated mainly in terms of annual change of cosmic ray diurnal 
variation. The results obtained are as follows: 

1) Annual change of cosmic ray diurnal variation is far more remarkable 
in the Antarctic region than in the moderate latitude. 

2) The best correlation between cosmic ray intensity and isobar level height 
is found at 50 mb in the Antarctic while at 300 mb in the moderate latitude. 

3) The amplitude of diurnal variation of 50 mb isobar level height at Syowa 
Base is estimated from the cosmic ray data at about 12 meters in summer and 


about 7 meters in winter, respectively. 
1. Introduction 


According to the results obtained from a Neher type ionization chamber data during 
the period from March 1957 to January 1958 at Syowa Base, Antarctic, it was found 
that the amplitude of cosmic ray diurnal variation was much larger in winter than 
in summer (Kitamura, Fukushima and Kodama, 1959). From the large seasonal change 
in diurnal variation of the surface temperature simultaneously observed, it is supposed 
that most of the above seasonal change in cosmic rays might be due to the atmospheric 
effect. On the other hand, the cubical meson intensity observed at Mawson during 
the same period did not reveal such a definite seasonal change as seen at Syowa Base. 
Therefore, one question remained concerning which of instrumental and meteorological 
conditions affects cosmic ray diurnal variation. 

In order to investigate the above question and also other related problems, a 
cubical meson telescope using plastic scintillators was constructed (Kawasaki and 
Miyazaki, 1958) and installed at Syowa Base on March 1959. The data obtained during 


March-December 1959 are analysed in this paper. 
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2. Annual Change of Cosmic Ray Diurnal Variation 


It was already found that the surface temperature at Syowa Base revealed larger 
diurnal variation in summer than in winter. This fact leads us to suppose a remarkable 
annual change in cosmic ray diurnal variation. To ascertain this supposition, the cos- 
mic ray data from three stations in the Antarctic region (Syowa Base, Mawson and 
Wikes) were investigated and two stations in Japan (Tokyo and Sapporo) for com- 
parisons. The first harmonics of diurnal variations (we describe hereafter “the amplitude 
of Ist harmonics of diurnal variation” as “ diurnal variation” in the sake of simplicity.) 
for each month were calculated on the meson intensities corrected for barometric pressure 
at these stations. Annual changes of the diurnal amplitudes at the Antarctic stations 
are given in Fig. 1 (A) and also those at Japanese stations in Fig. 1(B). It is clearly seen 
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Fig. 1 The annual changes of the diurnal amplitudes for the 
meson component. 

(A) Higher latitudes (Mawson, Wilkes and Syowa Base: 
Antarctic stations). 

(B) Moderate latitudes (Sapporo and Tokyo; Japan). © 


from the figures that the apparent annual change of the amplitude appears in the 
Antarctic region, while no changes appear in the moderate latitude such as Japan 
Generally speaking, the cosmic ray meson intensity on the ground should be affected 


| 
| 
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by both the primary cosmic rays and the upper atmospheric temperature. But the 
annual change of diurnal amplitude for the neutron component does not show any 
apparent difference between the higher and the moderate latitudes, as seen in Fig. 2. 
Therefore, the remarkable annual changes in Fig. 1 (A) seem to be not due to the 
instrumental but to the atmospheric effect. 


0-5 Mt.NORIKURA -. MAWSON BASE 


| je 34--ST6 7 S19 I10°T) 12 
MONTH 
Fig. 2. The annual changes of the diurnal amplitudes 
for the neutron component at the higher (Mawson 
Station, 1959) and moderate latitude (Mt. Norikura, 1959). 


3. Atmospheric Temperature Effect 


From the facts mentioned above, a considerable difference in the atmospheric 
condition can be expected between the Antarctic region and the moderate latitude. 
Fig. 3 shows a difference between the maximum temperature and the minimum through- 
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Fig. 3. The annual ranges of temperaure (given by 
height representation) throughout a year at various 
levels of isobar. 


out a year at a number of isobar levels. The maximum-minimum difference becomes 
constant above 50 mb isobar level at Syowa Base (Seino and Kawaguchi, 1961), contra- 
ry to 300 mb in the moderate latitude. This fact suggests that a true temperature 
coefficient for cosmic rays at Syowa Base can not be determined unless the isobar 
levels from 300 mb to 50 mb are taken into account. Numerical values of negative 
temperature coefficient calculated for a number of isobar levels are shown in Fig. 4. 
They are deduced from a triple correlation among I (daily mean value of meson in- 
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tensity at Syowa Base), H (height of isobar level at Syowa Base) and N (neutron in- 
tensity at Mawson), where N is adopted in place of the primary cosmic rays. In Fig. 4, 


SYOWA BASE 1959 
Tinew 


ye 


a 
of nce BASE 1957 


700 500 300 200 100 S50 30 mb 


Fig. 4. The negative temperature correlations (linked line, 
right side scale is available) and coefficients (full and dotted 
lines, left side scale is available) at various levels of isobar. 


the best correlation between I and H is found at 50 mb isobar level, where a tempera- 
ture coefficient 5 %/km is obtained, being the same as in Japan (Wada, 1951). 


4. Diurnal Variation of Isober Level inferred from Cosmic Rays 


Now let us estimate the diurnal variation of isobar level at Syowa Base from 
cosmic ray data. In order to extract the part due to the atmospheric temperature 
effect from the cosmic ray diurnal variation observed on the ground, the annual change 
shown in Fig. 1 (A) must be firstly corrected for the primary cosmic ray variation, using 
the data of neutron component at Mawson. In practice, the diurnal variation in the 
primary cosmic rays does not always appear in the same fashion for both the meson 
and the neutron components observed on the ground, because of the difference between 
their yield functions in the atmosphere. Taking this point into consideration, the 
approximation can be assumed: 


Ale late la ds Fovrw Dap) aie (1) 


n 


where I is the amplitude of the first harmonic of diurnal variation in cosmic rays and 
I is its annual mean, » and m represent neutron and meson components respectively, 
and 7 the month. 4I is the part due to the atmospheric origin. The calculation has 
been carried out for each of six months which are divided into two groups: winter 
season (March, April and May) and summer season (October, November and December). 
These results are given in Table 1. Using both the calculated 4I and a temperature 
coefficient 5%/km obtained in the preceding section, it is possible to estimate the 


amplitude of diurnal variation of isobar level. height, resulting about 13 meters in 
summer and about 7 meters in winter. 
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Table 1. Calculations of Equation (1) for each month 


month Ig ate 4; | dl,/5%/km 
nr 
Apr. 1.9 %o 2.5 %o —0.6 %o 12% 
May. 1.9 1.6 0.3 6 
Jun. 22 2.0 0.2 4 
mean ig 
Oct. 2.0 1.4 0.6 12 
Nov. i 1.4 0.3 6 
Dec. 0.8 Ve —0.9 20 
mean 13 


5. Discussions 


To justify the estimated amplitude of diurnal variation of isobar level, the radio 
sonde data obtained in the moderate latitude and in the north polar region were ex- 
amined. Monthly mean values of diurnal variations of isobar levels at Tatenot+ and 


Kap Tobint* are given in Fig. 5, where they are derived from the values of four 


WINTER SUMMER 
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109m 


Fig. 5. The observational amplitudes of diurnal varia- 
tion of isobar level at Tateno (100 mb, 300 mb), Kap 
Tobin (100 mb) and the calculated ones at 
Syowa Base (50 mb). 
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+ Tateno: 36-03N, 140-08E, Sea level, Japan. The data during 1959 were adopted. 
++ Kap Tobin: 70-25N, 21-58W, 42 meters above sea level, East coast of Greenland. The data 
during 1958 were adopted. 
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observations a day at Tateno and twice observations at Kap Tobin. 

As seen in Fig. 5, the seasonal change in diurnal variation of isobar level is not 
evident at Tateno, at both 300 mb and 100 mb, but clear at Kap Tobin. These meteoro- 
logical phenomena are consistent with the fact that the large annual change of diurnal 
variation in cosmic ray meson intensity is observed at Syowa Base whereas any change 
scarecely observed in Japan. Therefore, it would be reasonable to consider that the 
apparent annual change of cosmic rays in the Antarctic region is mainly due to the 
atmospheric temperature effect. 

However, another question related to this point is the result obtained at Mawson 
where the annual change mentioned above is not so remarkable as at Syowa Base (Fig. 
1(A)). Although a definite explanation for the discrepancy between the result at Syowa 
Base and that at Mawson have not been made yet, it is plausible to assume that it may 
be due to the difference in the meteorological states between both stations. In fact, some 
evidence supporting this consideration is seen in Fig. 2, but sufficient data to discuss 
this problem in detail have not been obtained yet. Detailed study will be given in 
future. 

In concluding, the authors wish to express their sincere thanks to Dr. Y. Miyazaki 
and prof. Y. Tamura for their kind guidances and encourogements throughout the 
work, and also to Dr. M. Wada for his valuable discussions. Thanks are due to the 
Japan Meteorological Agency who kindly provided the radio sonde data for us. 
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The Atmospheric Electric Field on the Summit of Mt. Hiei (860 m) 


By Akira OKAWATI 


Geophysical Institute, Kyoto University 
(Read May 16, 1960; Received June 20, 1961) 


Abstract 


On the summit of Mt. Hiei (860m) in Kyoto, as well as at the Geophysical 
Institute of Kyoto University (70m) the simultaneous and continuous measurements 
of atmospheric electric potential gradient were made in the period from December 
1955 to March 1956, and the characteristic of the diurnal variation of the atmos- 
pheric electric potential gradient on the summit of the mountain in winter season 
was studied. For the purpose of more precise investigation, the measurements 
were made again in the period from February 1959 to February 1960. Both the 
diurnal and the seasonal variations were clarified. On the summit of Mt. Hiei, 
even in winter, the diurnal variation was not similar to that over the ocean— 
unitary variation—, not like that on the Alps and Mt. Fuji, but it was affected 
by lower atmosphere in daytime. It was found that the mean value of the potential 
gradient on Mt. Hiei was greater in summer than in winter contrary to that at 
the Geophysical Institute. 


1. Introduction 


By the studies of many investigators, the characteristic of the calm day atmospheric 
electric field at various places on the globe have become clear. Over the ocean, the diur- 
nal variation of atmospheric electric potential gradient has the mode varying in the same 
phase all over the world, depending upon the universal time. But on the land, diurnal 
variation has the locality after the geographycal and meteorological condition of the 
observing station. It is generally accepted that the variation of atmospheric electric 
conductivity is little over the ocean and the variation of atmospheric electric potential 
gradient is mainly due to the electric potential of the equalizing layer which is maintained 
and varied according to the world-wide thunderstorm activity, and on the land the atmos- 
pheric electric conductivity and the columnar resistance vary with the concentration of 
air-pollution and the activity of austausch which depend on the geographical and 
meteorological condition of the observing station, and so the diurnal variation of atmos- 
pheric electric potential gradient has the mode depending upon the local time. Especially 
in town, where is much air-pollution, the atmospheric electric field has large mean 
value and broad amplitude of diurnal variation, and it has very local character. This 
was also mentioned by the author in his former paper (1961). 

H. Israél, H.W. Kasemir, and K. Wienert (1951) made measurements of the electric 
potential gradient and air-earth currents at Buchau, Zugspitze (2960 m) and Jungfraujoch 
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(3470 m). They found that in summer the relations between potential gradient, air-earth 
current and conductivity, and their variation with time of day are similar to those of land 
stations; but in autumn the similarity is to the phenomena of ocean. They 
payed attention to the role of austausch to the elements of atmospheric electricity. In 
summer, the austausch extends up to the height of the station and brings with it the dis- 
turbances due to air-pollution, but there are no effects of air-pollution from night to early 
morning. In this period of day the austausch seems to be lowered below the height of the 
station. The variation of atmospheric electricity at the summit of high mountain thus 
depends on the variation of exchange layer caused by the austausch. H. Israél measured 
the elements of atmospheric electricity at both the Jungfraujoch (3470m) and the Sonnblick 
(3100 m) about 400 km. apart from each other, and found that the diurnal variations at 
both stations are quite similar to each other. This fact means that both stations are 
free from the effect of austausch layer, whereas the plain land stations have usually 
their own types of diurnal variations when they are considerably distant from each 
other. The stations on the summit of high mountains have the diurnal variation of 
atmospheric electric potential gradient just like that over the ocean, as there is no 
cause of the variation of atmospheric electric conductivity and columnar resistance. 
R.C. Sagalyn and G.A. Faucher (1954) made measurements of atmospheric electric con- 
ductivity and large ion content of air by the aeroplane, and they found the number of 
large ions often to vary irregularly with height within the austausch region. Above 
the upper boundary of the austausch region, the number of large ions becomes very 
small. They made a number of aircraft measurements and found the depth of ex- 
change layer having much large ions and small electric conductivity to vary diurnally 
and seasonally from 300 metres up to 3000 metres. 

The effect of the variation of exchange layer upon the atmospheric electricity is 
very significant. From this point of view, the simultaneous and continuous measure- 
ments of atmospheric electric potential gradient were carried out at both the summit of 
Mt. Hiei and the Geophysical Institute of Kyoto University in Kyoto City. 


2. Results of Observation 


As for the method of observation, it was explained in the former paper by the author. 
The measurements were made in the first period between December 1955 and March 1956, 
then the plane reduction factor or the absolute value of potential gradient at Hiei was 
decided by the data used for the analysis of minor shower cloud electrification. And 
so the absolute value at Hiei in the first period has no significance, it connot be com- 
pared directly with the value in the second period or those of other stations. In the 
second period of observaton, the plane reduction factor at Hiei was decided by the 
comparison with the absolute measurements made in the motor-pool close to the ob- 
serving station with the area of about 3000 m2 (40mxX80m). But this value means only 
the value of a special place on the top and has a little significance. 

The mean value of the calm day diurnal variation of atmospheric electric potential 
gradient at both Hiei and the Geophysical Institute of Kyoto Uniuersity in the first 
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Fig. 1. a map showing the positions of observing stations and industrial and thickly 
populated area of hatching. 


period is shown in Fig. 2. The local component curves which are got by subtracting 
the unitary variation from the mean curves are also shown. The diurnal variation of 
atmospheric electric potential gradient at Hiei is very different from that at the Geo- 
physical Institute. The amplitude of diurnal variation is much less at Hiei than at 
Geophysical Institute. Comparing the local component curve at Hiei with the other, 
it is found that the curve at Hiei has the maximum towards 12”, but at the other has 
it at 8", and the time when the potential gradient begins to increase is at 8” at Hiei 
but it is at 5% at the other. The curve at the Geophysical Institute has the ninimum 
of diurnal variation at 4”, but the value at Hiei, to the contrary, is larger than before 
and after that time. On the other hand, the unitary variation has its maximum at 4” 
(G.M.T. 19”). These facts show that the atmospheric electric field at Hiei is much more 
controlled by the variation of the electric potential of equalizing layer than at Geophysi- 


cal Institute. In the period of day between 0” and 8", the diurnal variation at Hiei 


is very similar to the unitary variation. This shows that the summit of Mt. Hiei 


(860 m) is above the upper boundary of exchange layer in this period of day in winter, 
At Hiei the diurnal amplitude of the local component curve is less than at the Geo- 
physical Institute, the time when the potential gradient begins to increase is three 
hours later than at the Geophysical Institute and the time of the first maximum is two 
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or three hours later than that at the Geophysical Institute. From about the ait 
pollution begins to increase in Kyoto City (or on plain land) asl poets gradient 
also begins to increase at that time, but the summit of Mt. Hiei sag is about 800 
metres higher than Kyoto City, is still above the exchange layer at the time se: suffers 
no effect from pollution. The austausch becomes active in advance of the time of sar 
and soon the upper boundary of exchange layer becomes higher. The air-pollution is 
conveyed up to the summit about two hours later than the plain aang therefore the 
potential gradient at Hiei has the mode of diurnal variation beginning to increase towards 
8" and having its maximum at 12%. In the period of day between 8” and 12”, the 
potential gradient at Hiei has the diurnal variation considerably different from the 
unitary variation, but it has not so broad amplitude of diurnal variation as at the 
Geophysical Institute. The concentration of air-pollution is damped to be conveyed up 


to 800 metres height. 
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Fig. 2, The diurnal variations of atmospheric electric potential gradient and the local component 
curves in the first period of observation between Dec. 1955 and Mar. 1956. 

Fig. 3 shows the frequency distribution of the occurrence of the values of potential 
gradient at each time. Within the part between 0” and 7", the frequency distributions 
are nearly normal and their maxima are about 80v/m, but thereafter they disperse 
broader. A curve is drawn along the upper side of the swelling of frequency of each 
time, and another one along the lower side. The lower curve is very similar to the 
unitary variatian of atmospheric electric potential gradient with the exception 
of the part between 10” and 16”, while the upper one differs from the unitary except 
the part between O” and 6”. The upper curve is rather similar to that at the 
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Geophysical Institute which has two maxima in a day. This means that on the days 
when the austausch is less active and the upper boundary of exchange layer is lower, 
the air-pollution affects Hiei only in the period of day between 0 vad 16”, but on the 
days when the austausch is active, the potential gradient at Hiei has almost the same 
diurnal variation as at the Geophysical Institute. 


Potential Gradient 
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Fig. 3. The frequency distributions of occurrence of potential gradient at each time of Hiei. 


In order to get more precise knowledge of the atmospheric electric field at Hiei, 
measurements were made again in the second period of observation between February 
1959 and February 1960, for a year. Almost all the same results as in the first period 
were got in winter, that is to say, the diurnal variation at Hiei is similar to the unitary 
variation in the period of day between O” and 6”, but it differs after 7”, and the time 
of maximum is about 12”, though the time at the Geophysical Institute is towards 9”. 
In summer, even at Hiei, the time when the potential gradient begins to increase and the 
time of its maximum are the very same as at the Geophysical Institute. The austausch 
is very active in summer and the boundary of exchange layer is far higher than 860 
metres; the summit of Hiei is almost in the same meteorological condition as the 
Geophysical Institute. The air-pollution which is generated on plain land seems to be 
conveyed very quickly up to the summit. 

As for the seasonal variation of atmospheric electric field, the mean value at the Geo- 
physical Institute is greater in winter than in summer as usual, but at Hiei it is smaller 
in winter than in summer to the contrary. At the summit of Mt. Hiei in winter, as the 
austausch is less active and it is often above the exchange layer, the air-pollution is 
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Fig. 4. The diurnal variations of atmospheric 
electric potential gradient and the local com- 
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Fig. 6. The diurnal variations of atmos- 
pheric electric potential gradient and the 
local component curves in Aug. and Sept. 
1959, in the second period of observation. 
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Fig. 7. The diurnal variations of atmospheric 

electric potential gradient and the local com- 

ponent curves in Oct. and Nov. 1959, in the 
second period of observation. 
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Fig. 8. The diurnal variations of atmospheric 

electric potential gradient and the local com- 

ponent curves in Dec. 1959, and Jan. 1969, in 
the second period of observation. 
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few and the potential gradient is smaller. But in summer the austausch is active and 
the air-pollution is much more than in winter, so the potential gradient becomes larger. 


The air on Hiei seems to contain the same degree of air-pollution as on plain land in 
summer. 


3. Conclusion 


At the summit of Hiei (860 m), the mode of the diurnal variation of atmospheric 
electric potential gradient is quite different from that at the Geophysical Institute of 
Kyoto University in Kyoto City, and also it is not the mode like those at the Alps and 
Mt. Fuji even in winter. In winter season, the summit is above the exchange layer only 
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in the period of day between 0” and 6”, and sometimes it is in the exchange layer only 
in the period of day between 10” and 16", The air-pollution is conveyed up there from 
the lower plain land for two or three hours in winter, but in summer it is done very 
quickly. At Mt. Hiei there is the seasonal variation of potential gradient with its 
maximum in summer and the minimum in winter contrary to that at the observing 
stations on plain land. 
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Abstract 


y-hematites with titanium (called for the sake of convenience in this paper 
y-titanohematite) were extracted from natural rocks. And the chemical, crystallog- 
raphic and magnetic properties were investigated. The occurrence of this 
ferromagnetic mineral is also discussed and the contribution of this mineral to 


the natural remanent magnetism of the rocks is described. 


1. Introduction 


As an important ferromagnetic solid solution responsible for magnetism of rocks, 
present author’s attention has long been focused upon the so-called titanomagnetite 
having cubic crystallographic structure; Fe,;0,-Fe,TiO,. His attention became doubled 
when Kawai, Kume and the present author found an exsolution phenomenon of the’ 
solution and published its solvus diagram (1954 a), because the exsolution of titano- 
magnetite offered a key to the solution for the causation of the reverse remanent 
magnetism of igneous and sedimentary rocks (1954 b). 

On continuing and extending his study, the present author has found existence of 
a solution which possesses structurally similar to the titanomagnetite but chemically 
and magnetically dissimilar ferromagnetic characteristics in several altered rocks such 
as propylite, green tuff, rocks containing dehydrated lepidocrocite etc.. The results of 
the chemical analysis of these minerals showed that they all are in a highly oxidized 
state of the titanomagnetite with exceedingly large ferric ferrous ratio. This oxidized 
state can well be understood in terms of the “vacancy of iron atoms” which exists 
in the crystal. Such minerals are called here y state titanohematite (or simply y Ti- 
hematite). It was also found that the intensity and stability of remanent magnetism 
of rocks containing such minerals vary in very wide range. Some of the remanent 
magnetism of these minerals may be explained on the basis of the chemical remanent 
magnetization as have been proposed by Haigh (1958) and Nagata et al. (1958), and 
some on the basis of Thellier’s isothermal remanent magnetism. 

On the other hand it has been well known that in many rocks there exist another 
solid solution with chemical composition same as that of the above-mentioned solid 
solution but with rhombohedral crystallographic structure (@ state). This solid solu- 
tion has been found by Pouillard (1950) who synthesized them and also by Kume (1955) 
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and the present author who found one of the members of the solid solution in some 
metamorphic and several deep-seated rocks in Japan. The solid solution has become 
very famous since Nagata, Akimoto and Uyeda (1952, 1953) found self-reversal mecha- 
nism of the remanent magnetism in Haruna dacitic pumice. Akimoto’s rigorous in- 
vestigations (1954, 1955) have contributed much to the crystallographic properties and 
magnetic characters of this series of solution. Uyeda (1955, 1956, 1957) has also made 
extensive study to account for the causation of the reverse remanent magnetism from 
various points of view. 

Despite of such vast experimental results mentioned above, very little data are 
known about natural occurrences of the y- and a Ti-hematite and the relative com- 
parison of their magnetism. On the other hand, regarding different roles of the three 
solid solutions playing in the remanent magnetism of rocks, a comprehensive under- 
standing seems to be lacking and, therefore, is strongly required at present. Con- 
sequently in this paper, the description on the results the present author obtained and 


those the other authors did is so arranged that they may give some more light on the 
above-mentioned problems. 


2. Method of Measurements 


Method of measurements by which the present writer made his experiment are 


Table 1 (Part 1, Volcanic rocks). 


No. of ) Geological 


: Ip li .| Mean |Mean 
Locality Rock samples} age (c.g.s./g) | (c.g.s./g) Ip/li} Decli. | Incli. 
Sakurajima, Andesite 48 Present |3.7x10-2|4.1x10-3] 9.0 | N 6E | +48 
Kagoshima 

Asama, Nagano Andesite 42 ditto 4.6x10-3|1.5x10-3} 3.1 | N 8E | +40 
(Tenmei lava) 
Tomuro, Hornblende 12 | Pleistocene |7.6x10-8|1.8x10-8| 4.2 | N 6W | +57 
Ishikawa Andesite 
Shigisan, Mica-garnet Pleio- 7 -3 10-1 2.3 |N sw] +52 
Osaka Andesite 19 Pleistocene 2.110 9x 
Tojimbo, Two-pyroxene | 45 | Pleiocene | 3.2x10-8|2.8x10-?| 1.1 | N 20E | +52 
Fukui Andesite 
Ojima, Aphanitic i .1x1078| 2.1x10-?| 2.9 | N170W | —60 
Fukui Andesite 110 ditto Six 
Hamaji, Sanukitic 20 ditto |8.7x10-*/1.9x10-*] 4.6 | N118E | —70 
Fukui Andesite 
Matsushima, | Hypersthene 35 ditto 1.2x1072|1.1x10-8|10.9 | N174W | —56 
Fukui Andesite 
Yoshizaki, Basaltic 22 ditto |1.1x10-*|3.5x10-8| 0.31] N  8E | +54 
Fukui Andesite 
Ikeda, Welded tuff | 37 | Miocene |2.3x10-*|4.1x10-*| 0.56] N 18E | +57 
Fukui Andesite 
| ee Upper -4 -4| 0.22] N 40E | +50 
bingy Liparite 10 | Creteecous | 2-110" | 9.3x10 
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given in the first place. The descriptions are devided into three parts, of which the 
one concerns the natural remanent magnetism and the thermal demagnetization of 
rocks. The X-ray, chemical, microscopic and thermomagnetic analyses of the rocks 
forming ferromagnetic minerals are given in the second. While, high pressure experi- 
ments with respect to the 7 Ti-hematite are done in the last part. 

(1) The natural remanent magnetism (N.R.M.) of oriented samples of various kinds 
of rocks were measured by a sensitive astatic magnetometer. Upper and lower 
magets were very strictly adjusted in such that the vectors of the both magnets can be 
fixed antiparallel. The precise descriptions of this apparatus appear in Kumagai and 
Kawai’s report (1953). The specimen to be measured was placed in front of the 
lower magnet and measurements were made in the geomagnetic field of the laboratory. 
In so doing both the instantaneous magnetic moment (Ii) due to the geomagnetic field 
applied to the specimen and the permanent moment (Ip) fixed to the rocks can separately 
From the results the K6nigsberger’s Q value, that is, Ip/li was 
determined for each rock specimen. The data of rocks containing titanomagnetite as 


be determined. 


a ferromagnetic mineral summarized in the columns of Table 1 and Figs. 1 (a)-1(f). 


Table 1 (Part 2, Plutonic rocks). 


; No. of | Geological I li » dene i 
Sipealit R ig p i irection of 
wise sakes ock samples age (c.g.s./g) | (c.g.s./g) Ip/Ii polarization 
: at d 
TRovama’ | Granodiorite | 9 | Falae%olc | 4. 1x10-5|8.7x10-5| 0.47 | in lower 
: hemisphere 
ts lea Granite 10 | Cretaceous | 3.4x10-5 |4.5x10-5| 0.76 | N 6W, +43 
Kad ; of 23 
Fuku’ | Granodiorite | 12 ditto | 5.2x10-*|7.3x10-| 0.71 | NI6E, +52 
Tsuruga nearly at 
ae Granite 1 ditto /1.1x10-|9.2x10-8| 0.12 | Fandom in 
hemisphere 
Takefu, rer : 
Ste Granodiorite | 48 wie 7.3x10~*|2.1x10-*| 0.34 ditto 
Niu, 
Fukui Granophane 13 Cretaceous | 2.01074 | 6.4x10-4| 0.31 ditto 
Yuzawa, ; : 
Niigata Granite 20 Tertiary |1.6x10-4|1.8x10-*| 0.89 | N 9E, +50 
Wada, : ; 
Nagano ditto 21 ditto 8.0X10-5 | 2.2x10-*| 0.36 | N21W, +41 
Unazuki, : Pal ; nearly at 
Toyama | Augen gneiss | 28 | T*7/S%0N |1.3x10-*|4.6x10-+| 0.28 | Tandom in 
: i hemisphere 
unazu, cee: 
Gyfa. Granodiorite 13 Permian | 4.2x10-5|1.3x10-4| 0.32 nearly 
ha : horizontal 
Hodatsu, ‘ 
Ishikawa ditto 32 ? 3.5x10-4|7.6x10-*| 0.46 | dispersive 
Shirakawa, ‘ 
Gyfu Granite 27 | Cretaceous | 6.9 10-4 | 1.1x10-3| 0.63 ditto 
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Thermal demagnetization of N.R.M. was undertaken on several rocks by heating 
them in vacuum under non-magnetic field. The sample to be measured was surrounded 
by a non-magnetic electric furnace with non-inductively wound platinum resistance wire 
and was placed in front of the astatic magnetometer. The measurements of N.R.M. 
and T.R.M. at various temperature were made on both the heating and the cooling 
processes. As has been reported by Kawai, Kume and the present author (1954 b), 
this method of measurement is effective to find out each fraction of the total remanent 
magmetism, in case the ferromagnetic minerals are not single phased, for example, 
and they have different Curie points. It becomes especially important when original and 
mother titanomagnetite has exsolved into her children titanomagnetites with different 
Curie points. In such case there arises a question in what direction and with what in- 


N 


(a) (b) 


(d) 
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(e) (f) 
Fig. 1. Direction of the natural remanent magnetization of rocks containing titanomagnetites as 
ferromagnetic minerals (in Schmidt projection). 
full circle: lower hemisphere 
hole circle: upper hemisphere 
(a): Sakurajima, Kagoshima (b): Asama, Nagano 
(c): Tomuro, Ishikawa (d): Tojimbo, Fukui 
(e): Ojima, Fukui (f): Yoshizaki, Fukui 


tensity each exsolved phase contributes to the total remanent magnetism (J,.) of 
rocks. At temperature T,, just above the lowest Curie point of the constituent phases, 
the remanent magnetism due to this phase (J,) disappears and the total moment at that 
temperature (J7,,) is the resultant vector of the rest phases having Curie points all higher 


than T,,. Therefore by subtracting J,,; from Jj the remanent vector of the lowest phase 
can be determined. 


i 
Ji=Sto—Jres 


When similar procedures be carried out at every temperature corresponding to 
the Curie points of the co-existing phases, each remanent vector, J, J>,..., Jn could be 
determined separately. 


Several representative examples of this measurement are shown in Figs. 2(A)-(F) 
and Fig. 3. 

(2) In order to carry out thermomagnetic analysis an appratus was designed com- 
bining an electromagnet with quartz spring balance. Very small amount of rock 
samples, say 0.1mg, can be measured. And the Curie point of the ferromagnetic 
minerals contained in the very samples are determinable by heating and cooling it in 
a temperature range from —180°C up to 900°C. The sensitivity of the balance is 
sufficiently high to measure paramagnetism as well as ferromagnetism of nearly all of 
the rocks in nature. 


For the X-ray analyses, chiefly “Norelco” X-ray diffractometer with Fe and Co 
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Fig. 2. Representative examples of the thermal demagnetization curves and their analyzed fractio- 
nal vectors of exsolved phases. 
J): Saturation magnetization of rock specimen at room temperature (c.g.s./g). 
J,: Saturation magnetization of rock specimen at any temperature (c.g.s./g). 
J,: Intensity of remanent magnetization of rock specimen (c.g.s./g). 


x: Weight percentage of exsolved titanomagnetite. 
x y: Weight percentage of exsolved titanomagnetite which is occupied by 
FI} y magnetization (multiplied by a factor p). 


+: One half of the fractional vector J, of exsolved phase with the Curie point 0,. 
J: That of with antiparallel direction to the above one. 


target was adopted to determine crystallographic structure and parameters of the 
ferromagnetic minerals. The lattice parameters of each minerals contained in rocks 


196 


were analyzed by means of the powder pattern charts. 
ferromagnetic minerals were undertaken, and chiefly the amount 
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Chemical analyses of some 
s of Fe’, Fe?” and 


Ti” ions in them were determined. The results concerning y Ti-hematite are also 
shown in the Table 2. 


Table 2. Results showing chemical composition, lattice constant and intensity 


(N.R.M.) 


0.5 


ie) 


Fig. 3. 


of saturation magnetization 


Thermal demagnetization curves of rocks 
containing ;-titanohematite. 


Saturation} Lattice 


No. of Chemical composition Chem. comp. 
speci-| Locality Rock (in_wt. %) (in_mol. %) _ mag. 
men FeO | Fe,0,| TiO, | MnO | SiO, | FeO |Fe,O,) TiO, | (C-G.S./g) 
1 vag x Propylite | 6.03|71.83| 5.65! 0.57 | 7.39 |13.88|74.42|11.70| 58 
ditto, 
2 | weathered | ditto 4.78|74.81| 4.51| 0.48 | 6.26 |11.24]79.20| 9.56) 61 
sand in situ 
Asuwayama,| Green 
3 Pula pio 13.58 | 70.66 |12.51| 0.19 | 2.50 | 23.99|56.14/| 19.87 50 
Syakudani, Green 
4 Fukui | lapilli tuft | 20-63] 64.61 | 10.51| 0.70 | 6.87 |21.61|59.16 | 19.23 48 
Awatabe, | Welded |... bled ireed caste penetra ba 
5 | Anatabe, | Welded | 10.03|65.10| 8.64| 0.42 | 4.32 |30.44|54.81|14.75| 63 
6 | Silla, | Quartz |96 79/7179] 1.13| 0 | 
Cctavstea. Tis ecanite : 55 | 0.89 |44.98]53.32| 1.70 89 
M ae 
4 Wabdien Propylite | 9.23|/71.12| 4.06] 0.49 | 6.84 |19.10|72.49| 8.41 71 
8* ditto ditto 4,95|73.89) 5.27| 0.64 | 8.18 |10.64|78.01|11.35 46 
ae Contact = 
Kamaishi- : 
mine, Iwate| ,if0n | 31.41|68.40| 0.0 | 0.14 | 0.78 |47.96|52.04| 0.0 90 
j deposits 


* This specimen is same as No. 7 but was heated to 300°C for 4 hours. 


constant 


(A) 


y-Titanohematites in Nature and the Role they play in Rock-Magnetism 197 


For microscopic observation of the polished section of rock samples, polarized light 
_ was used in order to distinguish the above-mentioned three solid solutions, because 
each solution possesses different reflective power, and each phase becomes determinable 
_by naked eye under the cross or parallel Nicol observation. 

The microphotographs, when necessary, were taken. These are illustrated in 


Table 3. Results of high pressure experiments 


Number Pressure atm. Temperature °C Resulting phase 
1 26000 220+5 a 
2 31000 160-45 a 
3 23000 150-45 a 
4 15000 200-5 a 
5 10000 100+5 Fi 
6 10000 150+5 7 22%+a 
7 27000 90+5 r3l1%+a 
J 


1.0 


0) s 
200 400 600C 


Fig. 4. Thermomagnetic curves before and after high pres- 
sure treatment. 

full circle: Original specimen of 7-titanohematite with 
a slight endomorphic titanomagnetite (having 
Curie point of 540°C). 

hole circle: Transformation of ;- to @ titanohematite is 
nearly completed after treatment ; the curve 
is fully comparable to that of the titano- 
magnetite initially contained. 

hole circle with dot: Transformation of 7 to a is imcom- 
pletely seen under this condition. 
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Fig. 5. Stability region of 7- and a titanohematite in relating 
to pressure and temperature. Properties of the starting 7- 
titanohematite is found in Fig. 4 and Table 2-No. 1. 
©: completely transformed from ;- to a phase under the 
applied condition. 
@: remained completely with no effect as 7 phase under 
the applied condition. 
©: partially transformed from 7- to a phase under the 
applied condition. 


Plates 1-12. In general very high magnifications were required and oil immersion 
method was frequently employed. 

(3) High pressure experiments were added to the above-mentioned measurements 
in order to find the region in which y- and a@ Ti-hematite is stable in nature. The 
apparatus used for this experiment is an anvil type squeezer and the pressure intensifier 
designed by Kawai and his co-workers at Kyoto University. The description of the 
details appears in the publication by Kawai (1957). The maximum pressure attained 
by an improved equipment of the former is 50,000 atm. at temperature 600°C. 

The pressures impressed upon y Ti-hematite are ranging from 0 atm. to 40,000 
atm. and the temperature not higher than 300°C. The results of this experiment, as 
shown in Table 3 and Figs. 4 and 5, indicate that the Stability region of a Ti-hematite 


expands and that of the y Ti-hematite contracts as pressure under which the rock in 
nature has been placed increases. 


3. General Description of y-Titanohematite 


Before going into the description of the present author’s examination, it is necessary 
to show the data of y Ti-hematite so far obtained by many authors. It is Baudisch 


<< 
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and Welo (1925) who first have found remarkable ferromagnetism in synthetic ferric 
oxide (Fe,O;) and distinguished it from ordinary hematite. They called the former + 
hematite and the latter @ hematite. On the other hand similar iron ore mineral con- 
_ taining considerable amount of TiO, was reported by Wager (1927) from Bushveldt 
igneous complex in south Africa. His finding seems to be of importance since it indi- 
cates that TiO, is soluble in y hematite. 

First synthetic works of 7 hematite which contains TiO, were done by Chevallier 
and Girard (1950). Crystallographic and magnetic observations were made on a series 
of solid solution between Fe;O, and cubic ilmenite, although those minerals whose 
chemical composition are close to the latter have not been synthesized yet. Basta also 
(1953) reported his precise crystallographic study on the mineral having cubic structure 
from Bon Accord in Transvaal. Cell dimension of 7 Ti-hematite containing 19.45 % 
TiO, was found by him to be 8.3475+0.0001 A. 

Nicholls (1955), on summarizing the data of Chevallier and Girard and also on 
reviewing the reports of the other authors, pointed out the probable existence of 
another cubic solid solution whcse terminal components are y hematite and y ilmenite. 

It is also worth noticing that a solid solution ketween Fe,O, and y Fe,O, was 
synthesized by Hagg (1935) and found a continuous decrease of cell dimension from 
8.380 A in magnetite to 8.322A in 7 hematite, suggesting that the Vegard’s law is 
satisfied. 

The solid solutions proposed by many authors are shown in Fig. 6 of ternary dia- 


TiO, 


FeO Fe,0, t-Fe,0, 


Fig. 6. Probable binary system of titanomaghemite previously 
proposed by many authors in the FeO-Fe,0,-TiO, ternary dia- 
gram. 
AB’: solid solution by Chevallier et al. 
BB’: solid solution by Nicholls. 
ABCD: 7-titanohematite having nearly trapezoidal com- 
positional field assumed by the present writer. 
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gram. On the diagram the thick lines BB’ represent chemical composition of Nicholls’ 
solid solutions and AB’ the Chevallier’s solution respectively. While the dotted line 
represents that of the titanomagnetite. The data of 7 Ti-hematite are so imcomplete 
and the descriptions are so fragmental that the present writer would like to make a 
comprehensive understanding with respect to the crystallographic, chemical, and magnetic 


behaviour in the following section. 


4, 7-Titanohematite as an Oxidized Product of Titanomagnetite 


As far as my results concern, the chemical compositions of all of the specimens are 
found in a triangular area ABC on the diagram of Fig. 6, and either cell dimension or 
intensity of saturation magnetization J, seems to vary rather continuously from place 
to place in this area. The changes therefores, are shown by equi-cell dimension and 
equal J, curves in Fig. 7, and also shown by numerical expression in Table 2. 
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Fig. 7. Relations between chemical composition and each of equal lattice para- 


meter and equal saturation moment at room temperature which are estimated 
from the present author’s specimens, 


From the above continuities and also from the present author’s consideration stated 
nereattets it seems highly plausible that y Ti-hematite is not always restricted to the 
binary system of solid solution as Nicholls and Chevallier et al. have proposed. The 
reason is that the mineral in question is generally found in nature as an oxidized 
and exsomorphic product of titanomagnetite as is revealed in microscopic observations 
of polished sections, so that its magnetic and crystallographic structure can more 
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naturally be explained on the basis of titanomagnetite possessing vacancies correspond- 
ing to the grade of oxidation. If the existence of the vacancies and its continuous 
variation from one to the other be the characteristics of the very mineral, then the 
chemical criterion in which the mineral possesses y-Sstate crystallographic structure 
should be illustrated by an area and not by the lines in the ternary diagram. 

Because magnetism of the mineral depends upon the ordering of the vacancies when 
they were many in number, one must consider an oxidizing process of titanomagnetite 
from which magnetic properties and crystallographic structure of the oxidized mineral 
can be inferred. 

In general, oxidation of ferrite always takes place in such that the Fe!” ions 
inside the ferrite are pulled out and are combined at the surface absorbing near-by 
existing oxygen ions with the instantaneous oxidation of the ion into Fe””. As a 
result there remained in the crystal cation vacancies whose number corresponds to 
that of the expelled Fe’’ ions. Electric neutrality is satisfied by changing two Fe! 
which are lying close to the above-mentioned vacancy into two Fe” ions. The situa- 
tion can be shown by the following. expression. 


(oxidation) 
3 Fe —— > 2Fe™ + [] + |e kal ear: 
remained vacancy expelled 
(crystal surface)  (“Smames) — (VRaa") (oxidized) 

1 

1 +44 ea an 
‘ oe 
a LEE a aha \ +t YO +tt 
S Sa Ss 


(inside of the crystal) 

In case of the oxidation of pure magnetite into y Fe,O, three Fe’? ions in the 
16c sites can be changed into two Fe” when one of them escapes out of the crystal. 
Therefore, in a unit cell, 8/3 times of such escape render the original magnetite to 
change into y-hematite leaving 8/3 vacancies therein. 

It is quite reasonable to assume that these vacancies make the cell dimension of 
y hematite to be slightly smaller than that of the original magnetite. Cell dimension 
a=8.39 A in magnetite decreases to 8.32 A when the former are changed into 7 hematite. 
Similar kind of oxidation would be valid even when the original magnetite contains 
Ti atoms in 16c sites of spinel structure. This corresponds to the case of the oxidation 


of titanomagnetite whose Ti content is less. 
5. Magnetism of ;7-Titanohematite 


Curie point measurement is very difficult to carry out with the present author’s 
specimens because during thermomagnetic analysis 7 Ti-hematite is transformed into 
a phase at temperatures ranging from 210°C to 500°C. Some representative thermo- 
magnetic curves showing their transformations into « phases are illustrated in Figs. 8 
and 9, and also, the dependency of transformation points of 7 titanohematites to their 
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chemical composition are shown in Figs. 11 and 12 respectively. 

While, changes of saturation moment y with the time of heating x of both the 
synthetic y hematite (A) and the natural 7 Ti-hematite (B) at their transformation 
temperature are illustrated in Fig. 10, and the rate of transformation 1—y at any time 
x can be given in the following formulae. 

1—y=1.01 e~%5*—0.03 (for A) 

1—y=0.92 e-0-110x""+ 0.14 (for B) 
The velocity of transformation is found by the differentiation of the formulae mentioned 
above. 

Regarding transformation from y- to a hematite as well as the transformation in 
y Ti-hematite, neither the detailed description nor the decisive conclusion has ever been 
seen yet. Several authors have only reported that the temperature at which the trans- 
formation takes place in Ti free y hematites is found to be ranging from 200°C to 
700°C. One of the answers in the present author’s mind is that 7 hematite is unstable 
in a wide temperature range (even near room temperature) but becomes stable when 
in the vacant seats of the mineral are embeded some foreign atoms, as the results 
rendering the very mineral to be stable even at moderately high temperture. 

The Curie point determination can be made only when the transformation tem- 
perature is raised higher than the Curie point of the mineral. When the above- 


Normalized Scale 


200 400 6006 
Fig. 8 Abrupt chages in saturation moment depending upon 
transformation at various temperature are found in some chemi- 
cally different kinds of specimens. 
full circles and line: heating process 
crosses and dotted line: cooling process 


@: magnetite crystal, Kamaishi-Mine @: Asuwayama, Fukui 

(2): Sdjiya, Okayama (6): Awatabe, Fukui 

(3: Syakudani, Fukui (6): Maruyama, Fukui 
~ , 


Fig. 10. The similar mode of progress in transformation of both synthetic 7 hematite 
and natural ;7-titanohematite can be seen from the time variation in each saturation 


A: y hematite synthesized by Takada (1960) kept at the transformation tempera- 
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Fig. 9. A thermomagnetic analysis confirming each y-titanohematite 
component with different transformation point and its endomorphic 
titanomagnetite. 
sample: specimen No. 1 in Table 2. 
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ture of 430-440°C°. 


y-titanohematite same as the specimen No. 6 in Fig. 8 kept at the transforma- 


tion temperature of 230-240°C°. 
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Fig. 11. Dependence of the transformation point of ;-titano- 
hematite either on its content of TiO, or on its numbers of 
metalic atoms to 32 oxygen; (Fe+Ti/O) x32. 

full circle and solid line: for TiO,-axis : 
open circle and broken line: for (Fe+Ti/O)x32 axis 


TT.FO, 


(Ti +70 x32) 


F, 0, 09(F /R+ ty? 


TRO, 


Fig. 12. Relation between chemical composition and transformation point of y titanohematite, 
represented on Fe/(Fe+Ti) vs. (Fe+Ti/O)x32 diagram. 


mentioned vacant seats are embeded by non-magnetic ions for example by Na/*, the 
transformation is suppressed and the Curie point becomes estimable. However, such 
data appear so few that further studies are required. Thus, the sure identification of 


* Michel and Chaudron (1935) determined Curie point of such 7 hematite and reported it to be 
675°C; 
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the mineral by means of thermomagnetic analysis is impossible. The only approach to 
the clarification of the magnetic character as well as the identification is the measure- 
ment of saturation magnetic moment J, by referring the crystallographic structure. 

Intensity of saturation magnetic moment of Ti free 7 hematite is expected to be 
5/2 #4, per mol. when we assume the expelling model of Fe’! to take place from pure 
magnetite as shown in eqation (1). This value, 5/2 “, per mol. in 7 hematite shows 
a fair agreement to the Chevallier’s model (1955) of 7 hematite. 

The question to be considered is what kind of structure and what magnetic moment 
these escaped Fe”’ atoms possess at the surface. 

(A) When they are formed into a hematite with absorbing oxygens surrounding 
the crystal surface, they should possess non-magnetic moment due to its antiferro- 
magnetic character. Therefore the bulk magnetic moment decreases from 5/2 4, slightly 
according as the surface layer of a phase becomes thick. 

_ But in case (B) the escaped ions are formed into y-hematite at the surface, the in- 
tensity suffers no change from the theoretical one. 

The change of magnetic moment due to the oxidation can more properly and more 
generally be expressed as follows; when the spin momenta of the expelled Fe” ions 
would supperimposed upon those of the Fe” ions inside of the crystal, then the increase 
of the magnetic moment could be expected (C). On the other hand, if the former 
momenta would off-set the latter, large diminution of magnetic moment should occur 


2.0 


LO 


0 10 20 30 40% 
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Fig. 13. Expected relation between saturation moment of fully 

oxidized ;-titanohematite and its content of assumed 7TiFeO, 

(y). 

A: calculated values of completely oxidized 7-titanohematite 
from Akimoto’s model of titanomagnetite through the present 
writer’s expelling mechanism. 

B: that from Chevallier’s model of titanomagnetite through the 
present writer’s expelling mechanism. 

C: calculated saturation moment of oxidation product from 
original titanomagnetite after S. Akimoto (1957). 

solid circle: 7-titanohematites specimens nearest to the system 

of 7F,0;-7TiFeO, (refer to Fig. 7). 
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(D). Details of the matter, however, should be discussed in term of “ superexchange ” 
of 3delectrons through the oxygen in metalic oxide, especially, in the neighbourhood 
of the boundary layer of the original and oxidized titanomagnetite, as is often dis- 
cussed in very fine powder or thin film. The present author’s results demonstrate that 
(A)-type oxidation took place and the intensity of saturation moment, and cell dimension 
have always been decreased in the case of oxidation of titanomagnetite which contains 
rather small amount of TiO, (Fig. 13). 

The above-mentioned vacancy production and magnetic change which occurs in 
16c sites, however, are inadequate in case of titanomagnetite whose Ti concentration is 
very high. In such a mineral Fe” ions are not only placed in 16c sites but also in 
8f sites of which some part should become vacant as the oxidation goes on and 
on. To this extremity no exact interpretation has ever been given yet. Akimoto, 
Katsura and Yoshida recently (1957) obtained oxidized titanomagnetites which possess 
strong ferrimagnetism. So that, one must consider such an expelling mechanism of 
Fe” in sites of not only 16c but 8f sub-lattice that the original magnetic moment 


may increase. 


6. Remanent Magnetism and its Stability of 7-Titanohematite 


Since y Ti-hematite is the product of low temperature oxidation of titanomagnetite, 
or otherwise the product of the dehydration of lepidocrocite at moderate temperature, 
N.R.M. might be attributed to C.R.M. or to I.R.M. and not in general to T.R.M.. The 
reason is that the temperature at which the remanent magnetism has been acquired 
is lower than the Curie point of the mineral. 

So far as the present author’s examinations concern, N.R.M. of 7 Ti-hematite 
is generally unstable as compared to that of the titanomagnetite. It must also be 


Fig. 1, Change in the direction of the remanent magne- 
tization at any duration of time in the laboratory ; its 
causation of N.R.M. may be assigned to I.R.M.. 
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Fig. 15. Relation between directional change and 
duration of time (log. scale). No. of specimen corres- 
ponds to that of Fig. 14. 


Fig. 16. Natural remanent magnetism with good 

stability in Schmidt’s projection. 

hole circle: upper hemisphere. 

full circle: lower hemisphere. 
emphasized that the magnetic stability of 7 Ti-hematite itself varies in a wide range; 
some rocks have rather high stability as well as coercivity (Hc) in geomagnetic field 
and some have very low. Schmidt’s projections Fig. 14 show the change of remanent 
magnetism of rocks which are placed in the geomagnetic field of laboratory for 
sufficiently long duration. Whereas, in Fig. 15 the time change of direction of N.R.M. 
is shown. 

This method of judging the stability, strictly speaking, is not a quantitative one 
at all but is indeed a very practical way by which it is inferred qualitatively. Kawai 
and Kume (1953, 1959) assumed that remanent magnetism was unstable when they got 
a large change of direction and intensity of N.R.M. during the prolonged preservation 
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N 
5 
Fig. 17. Natural remanent magnetism with Fig. 18. Natural remanent magnetism without 
metastability. stability. 


with fixed orientation in laboratory. According to this judgement, the present writer 
classified the remanent magnetism of rocks containing y Ti-hematite into the following 
three parts ; (1) those with extremely high stability (Fig. 16), (2) those with metastability 
(Fig. 17) and (3) those without stability (Fig. 18). 

Reflective microscopic observations done for the three types of rocks demonstrate 
that the rocks having unstability of their remanent magnetism possess 7 Ti-hematite 
whose grain size is fairly large, say 0.01 mm in diameter (or width), and also that the 
rocks showing good stability of the remanent magnetism possess y Ti-hematite whose 
grain size is very small, say, several microns (see plates 6 and 7). The stability of 
remanent magnetism, especially its dependency upon grain size can quite reasonably 
be explained either on the basis of Néel’s thermal fluctuation theory for single domain 
particles (1949) or on the basis of the theory of domain wall displacement (1951 a, 1951 
b, 1955). 

Nevertheless, a question still remains as for the afore-stated fact that + Ti-hematite 
is found to possess less stability than the stability which the corresponding titano- 
magnetite having same grain size possesses. One of the most probable answer might 
be given by the “existence of vacancies” in y Ti-hematite. Nakamura (1957) recently 
found in metalic iron an effect of vacancy upon the domain wall movement. The 
specimen once heated at high temperature and quenched, and therefore, possessing 
large numbers of vacancy at room temperature showed remarkable magnetic after-effect. 
He considered that the vacancies enable the domain wall displacement to occur easily. 
If such a kind of interpretation be available also in oxides such as y Ti-hematite, the 
above answer could be accepted with less resistances. 

Generally the vectors of N.R.M. of the rocks containing this mineral produced by 
hydrothermal reaction is found frequently to be at random in both its direction and 
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the intensity. Even when one of the rock samples is cut into several pieces, it is 
observed that each piece possesses different vector of remanent magnetism. This seems 
to be chiefly due to the heterogeneous distribution in the very rocks of stable and 
unstable remanent magnetisms. This tendency is rather characteristic when the altered 
titanomagnetites with heterogenious distribution in rock possess unusually strong intensity 
of remanent magnetism (Fig. 19). When such rock is magnetized in a uniform 
magnetic field, static magnetic energy would become increased so enormously that one 
directional magnetization becomes unlikely. We can expect such a flux closure lines 
of force which decrease the static energy. 


Fig. 19. An unusual distribution of direction and intensity in N.R.M. im situ. 
values beside the circle mean intensity with the order of 10-? 
C.G.S./g.; (+) values and (—) values indicate the inclination 

to downward and upward respectively. 


N 


s 
Fig. 20. Remanent polarization of rocks containing original titanomagnetites having been 
free from hydrothermal alteration is illustrated referring to that of the Iwaine 
Formation, which has the same geological age of middle Miocene. 
@: Itd Formation, Fukui Prefecture (by the present author). 
x: Iwaine Formation, Toyama Pref. (by Kawai). 
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On the other hand, even in the same series of rocks where no hydrothermal re- 
action, hence, no alteration of titanomagnetites has occurred, the direction of the rema- 
nent magnetism is found to be uniform everywhere in the specimen and at the same 
time their magnetism is fairly stable. This kind of magnetization can be considered 
in most cases to be T.R.M.. Age of these rocks are assumed from geological evidences 
to be middle Miocene. In 1950 Kawai (unpublished): observed N.R.M. of rocks of 
similar age from Iwaine formation, in the Central Japan. It seems to be worth notic- 
ing that his results show a fair agreement to those by the present author (Fig. 20). 


7. Stability of 7-Titanohematite 


The significance of carrying out high pressure experiment with 7 Ti-hematite lies 
in the following facts. In some metamorphic or plutonic rocks which were formed 
under reasonably high pressure, one can rarely find y Ti-hematites. Nicholls* 
predicted the existence of cubic ilmenite in some rocks which have crystalized under 
very high pressure. This is a remarkable contrast to the present author’s observation. 
Under such high pressure y Ti-hematite, even if it had formed, should have changed 
into @ Ti-hematite. It is evident that the relative volume of y Ti-hematite is about 
7% larger than that of a Ti-hematite, so that the decrease of Gibb’s free energy 
due to PV term becomes greater when the former is transformed into the latter (a 
phase) under high pressure. This tendency is especially remarkable in some metamor- 
phic and deep-seated plutonic rocks. 

According to the recent synthetic work on this mineral, y Fe,O, is metastable 
substance. It exists only in a transient chemical reaction such as oxidation of titano- 
magnetite or the dehydration of y FeOOH. The life time of the mineral depends 
upon both the grain size and the amount of the impurities and becomes longer as the 
latter increases. 

As the present author’s experiments indicate (Fig. 5), under both the high pressure 
and the temperature the region in which the minerals exist becomes narrow. 


8. Natural Occurrence of 7-Titanohematite 


yr Ti-hematite is by no means a generally occurring mineral as compared to 
titanomagnetite and «@ Ti-hematite. The latter two solid solutions are very common 
in volcanic, plutonic and pyroclastic sediments as have been reported by nearly 
all of the authors whose subjects of the study concern rock magnetism. Whereas, the 
rocks containing 7 Ti-hematite are limited to some igneous rocks such as propylite, 
altered green tuff and some sediments. As those petrological properties indicate, the 
rocks have more or less been altered after their original formation by the subsequent 
hydrothermai reactions. 

The original titanomagnetite might have been oxidized directly into y Ti-hematite 
perhaps with low temperature and probably with sluggish reaction. This is a kind of 


very slow rate oxidation of the rock or soil containing ferromagnetic minerals during 
their weathering. 
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Another point to be emphasized is that y Ti-hematite can also be produced by a 
dehydration of lepidocrocites (7 FeOOH) when they are heated at moderate tempera- 
ture ca. 200-400°C. This kind of formation have been confirmed by the pressent author 
during his thermomagnetic analyses. On heating the rocks containing lepidocrocites 
at temperature ranging from 150°C to 250°C, ferromagnetism appears with or without 
simultaneous dehydration of the mineral. On further heating up to 550°C the ferro- 
magnetism once appeared disappears. The appearance of ferromagnetism is indeed 
due to the production of ; hematite and the disappearance to the decomposition of the 
ry phase into a phase. 

A method of producing 7 hematite by means of synthesized lepidocrocite has 
recently been found by Takada (1960) who also found the change of 7 FeOOH into a 
Fe,O; through the intermediate product of 7 Fe,O3. 

In summarizing the manner in which y Ti-hematite is produced, following three 
cases are of importance in nature: 

(1) Low temperature oxidation of titanomagnetite by hydrothermal reaction. 

(2) Low temperature and slow rate oxidation of titanomagnetite by weathering. 

(3) Dehydration of lepidocrocite. 

As has been already mentiond 7 Ti-hematite is impossible to be found from deep- 
seated and metamorphic rocks unless they have suffered secondary alteration. This 
fact is supported by the data of the high pressure experiments. Consequently, we can 
conclude that the occurrence of 7 Ti-hematite is rather limited to the pressure free 
and oxygen rich part of the earth crust, for example, in the uppermost surface layer 


P of the crust. 
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7-Titanohematite completely oxidized from titanomagnetite, (x 400). 

Light grey: 7-titanohematite. 

Grey: slightly remained titanomagnetite. 

Black: pits and cracks. 

Rock: green tuff, Asuwa, Fukui Pref.. 

Enlarged pattern of plate 1, (x1200). 

Gradual oxidation of titanomagnetite to 7-titanohematite, (800). 

White: a-hematite. Light grey: jy-titanohematite. 

Grey: titanomagnetite. Rock: welded tuff, Awatabe, Fukui. 

7-Titanohematite nearly oxidized from titanomagnetite, (1000). 

Rock: green tuff, Asuwa, Fukui. 

Remark: oxidation of titanomagnetite to ;7-titanohematite commences at cracks or 
partings of the former and proceed towards the inside. 


Another photograph of the same rocks as plate 4. Regular development of exsomorphic 
y-titanohematites with inherent direction are noticeable, (x 1000). 

Enlarged photograph of the same sample as plate 5, (2000). 

Minute veinlet precipitates of 7-titanohematite, (x 2000). 

Rock: propylite, Maruyama, Fukui. 

Remark: the rocks possess stable N.R.M.. 

Another pattern of 7-titanohematite in titanomagnetite, (600). 

Rock: green tuff, Syakudani, Fukui. 

Remark: the rocks possess unstable N.R.M.. 

Original structure in titanomagnetite (phenocryst) and the characteristic partings therein, 
(x 1200). 

Rock: propylite, Maruyama, Fukui. 

Photograph showing the progress of oxidation along minute partings, (1800). 

Rock: propylite, Maruyama, Fukui. 

Exsolution lamellae of ilmenite (relict) in titanomagnetite phenocryst, (400). 

Rock: propylite, Maruyama, Fukui. 

Local and intensive oxidation directly into a-hematite can be seen, and intermediate 
product of 7-titanohematite can only slightly be observed, (600). r 

Rock: propylite, Maruyama, Fukui. 

White: a-hematite. 

Grey: titanomagnetite. 
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Abstract 


Cosmic ray neutron dada obtained at Syowa Base, Antarctic, are analyzed for 
studies on characteristics of the four cosmic ray unusual increases, which occurred 
on May 4, 1960 and during November 1960, by comparing with cosmic ray data 
from the other stations in the polar regions. Investigations of the enhanced 
intensity, the onset time and the decay time of the cosmic ray increase give 
definite differences between the northern and the southern polar regions, differing 
from event to event. Such differences suggest that the polar type increase is 
not only raised by the solar cosmic ray particles distributed isotropically in the 
interplanetary space but also by the anisotropic particles uniformly deflected from 


their initial emitted direction due to the solar magnetic cloud. 


1. Introduction 


Generally speaking, solar cosmic ray particles ejected from the sun reach the earth 
in manner of either direct or indirect propagation through the interplanetary space. 
The direct incident particles, which persist their initial emitted direction to the 
vicinity of the geomagnetic field, give rise to the so-called impact zones on the earth. 
The indirect particles appear as increase event in the polar regions above 70° geomagne- 
tic latitude where the incident are usually forbidden for the direct particles (Jory, 1956; 
Liist, 1957). In the latter case, if solar cosmic ray particles become completely isotropic in 
the space, a similar increase would be possibly expected at every station in the polar 
regions. However, if not so, it can be supposed that there may be some differences 
in characteristics of increases among these stations. Such differences would also 
arise by assuming such a modulation effect on the direct component as deflection due 
to the magnetic cloud emitted from the sun (Pomerantz et al, 1961). 

Recently, unusual increase events of cosmic rays frequently occurred, that is, once 
on May 4, 1960, twice on November 12, 1960 and once on November 15, 1960. For- 


* S. Fukushima was lost in a fierce blizzard in the Antarctic on October 17, 1960. This is one of 
his posthumous works. 
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tunately, as several cosmic ray neutron monitors had been operating at the southern 
polar stations during 1960, the direct comparison of their data with those from the 
northern polar stations is possible for each event. The neutron monitor installed 
at Syowa Base (69.0°S, geomagnetic) on March 1960, too, recorded these increase 
events which data (Appendix) are available under the above consideration, especially, 
they are analyzed together with the data from Churchill (68.7°N, geomagnetic) corres- 
ponding to the opposite latitude of Syowa Base. Detailed discussions using the data 


from other cosmic ray stations in the both polar regions are given. 


2. Observational Results 


The type of neutron monitor at Syowa Base is not identical with the standard type, 


but cosmic ray neutrons measured by it show time variations similar to those by the 


standard type (Fukushima, Kitamura and Kodama, 1961). For direct comparison 


of the data from Syowa Base with those from Churchill, both 15-minute values 
of the former and 10-minute values of the latter were reduced to the relative intensities 
by taking the pre-increase intensities as the normal intensity levels. Time variations 


of the relative values for four events are illustrated in four diagrams of Fig. 1, where 


" 12 «18 
MAY 4, 1960 NOVEMBER 15, 1960 


Fig. 1. Comparison of neutron intensities among two stations in the four unusual 

increase events of cosmic rays. Rigid and dotted lines are for Syowa Base and 

churchill, respectively, where the dotted line before the double vertical line in the 
November 12th event is for data from Ottawa. 


rigid and dotted lines are for Syowa Base and Churchill, respectively. In the case of 
the first event of November 12, data from Ottawa are adopted in place of data from 


Churchill which were in failure at that time. 


Prior to consideration about characteristics of these increase events, let us first 
examine the differences between two stations, Syowa Base and Churchill, concerning 
the following properties: the maximum enhanced intensity, the onset time of the in- 


crease, the rise time and the decay time of the enhanced intensity. The rise time is the 
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difference between the starting time of the intensity enhancement and the time of the 
maximum intensity. The decay time means the interval during which the enhanced 
intensity decreases by half. Comparisons of the above properties at Syowa Base related 


to those at Churchill are. summarized in Table 1. 


Table 1. 
1) May 4 | (2) Nov. 12| (3) Nov. 12| (4) Nov. 15 
Event Mi Ist 2nd 
Ratio of Intensity 5 1 1 1 
Difference of Onset Time 0 0 _ +1 hr 
Ratio of Rise Time 1 1 = 0.5 
Ratio of Decay Time 1 2 1 1 


Next, another comparison was made among the all stations belonging to the both 
polar regions: Thule, Resolute, McMurdo, Mawson, Ellsworth, Syowa Base and Chur- 
chill. In that case, the maximum enhanced intensities were deduced from 15-minute 
values for every station except Churchill where 10-minute values were reported, and 
then, for convenience, normalized to the value from Thule. Ratioes thus obtained for 
four events are given in Fig. 2, resulting a large fluctuation in the enhanced intensities 
in the case of Event 1 alone. Comparing these results from event to event about the 
rise time and the decay time, their visual properties common to all stations are 
summarized in Table 2, together with the results of Fig. 2. Physical meaning deduced 
from Tables 1 and 2 will be discussed in the following section. 


Table 2. 
No. of Event (1) (2) (3) (4) 
Rise Time shorter long | short short 
Decay Time short longer long long 
aon epi 24 0.5 0.2 0.5 


* Maximum range of ratioes of each station to Thule 


3. Discussions 


Most of solar cosmic ray particles starting from the vicinity of the sun can not 
reach the polar region above about 70° geomagnetic latitude on the earth’s surface, 
but if could, the intensity enhancement of cosmic rays observed there would be few. 

But in the case when the sun is situated in high latitude referring to the geomagnet- 
ic coordinate system, some polar regions belong to the impact zones. At the present 
events, where the source latitudes are 16°N for the May event and 19°S for the Novem- 
ber events, the stations shown by open circles in Fig. 2 belong to the impact zones 
while black circles for the polar regions, assuming a magnitude of the source as +15°. 

As the increase in the polar regions is usually caused by the isotropic solar cosmic 


ray particles, it is natural that apparent characteristics of the polar type increase should 


be nearly the same everywhere in the polar regions. In other words, numerical values 
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80° 
MAY 4th EVENT 


NOVEMBER 12, Ist EVENT 


NOVEMBER i2, 2nd EVENT 


RATIO OF MAXIMUM ENHANCED INTENSITY AT EACH STATION TO THAT AT THULE 


R Me Ma SC E 
Fig. 2. Comparison of the maximum enhanced intensities of the 
increases among seven stations belonging to the polar regions, 
where each station is arranged in order of the 
conventional geomagnetic latitude. 


about the various characters given in Tables 1 and 2 should be normally close to 1 
for the ratio and 0 for the difference and the intensity fluctuation. However, in the 


- case of the increase events occurred during 1960, the practical values about the above 


characters and the qualitative properties suggest that a similar type of increase does 
not always appear in the both polar regions and moreover characteristics of the in- 


’ crease differ from event to event. Phenomenological features of each event are as follows.. 


Event 1 (May 4th) 
Two stations, Churchill and Ellsworth, belong to the impact zones. It explains the 


enhanced intensities at the both stations are larger than the others, but the same rise 
time in everywhere is rather unreasonable. The most important feature of this event 
is that the decay time is shortest (Anderson et al, 1960). The fluctuation of the enhanced 
intensities still remains among the stations except Churchill and Ellsworth and at least 


it is larger than in the other events. 


_ Event 2. (November 12th, 1st) 


Mawson and Syowa Base are inside the impact zones, nevertheless, they do not 


display any features different from the others. A noticeable fact is seen about the decay 


time in Table 1. 
Event 3 (November 12th, 2nd) 
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No station belongs to the impact zones, being consistent with the fact that few 
fluctuation in all characters of the increase is found among all stations. According 
to the other analysis about this event, the apparent increase observed on the ground 
was caused by both the effect of distorted geomagnetic field and the part of solar 
cosmic ray particles trapped inside the solar magnetic cloud (Kodama and Kitamura, 
1961). Properties shown in Tables 1 and 2 do not contradict with such assumptions. 
Event 4 (November 15th) 

The onset time does not coincide among the both polar stations, i.e., about one 
hour earlier at Syowa Base than at Churchill. Similarly, the rise time at the former 
station is about one half at the latter. Ellsworth station belongs to the impact zones 
but displays no intensity increase distinguishable from those at other stations. This event 
happened during the development of a large Forbush decrease. The above time 
difference would be contributed to modulation effect of solar cosmic ray particles by 
the solar magnetic cloud by which the decrease was caused. 

It was pointed out by McCracken and Palmeira (1960) that the ratio of the rise 
time to the decay time is constant in all increase events, by examining on cosmic 
ray data of the background zone. But considering the polar type increase, it is deduced 
from Table 1 that the ratio suggested by them is not constant through Events 1, 2 
and 4 This fact suggest that the polar type increase, or solar cosmic ray particles 
entered the polar regions, contains the particles different from those entered the impact 
zones and the background zone. 

Now we suppose that the polar type increase consists of two kinds of solar cosmic 
ray particles. The one is the particles which become completely isotropic and uniform- 
ly fulfilled in the interplanetary space. The other is those deflected from their initial 
emitted direction but not still fulfilled in the space. If the polar type increase 
could be caused by the former alone, the north-south asymmetries as observed in the 
present events would not be expected. The particles permitted to enter directly below 
about 70° geomagnetic latitude will be deflected toward the polar regions, assuming the 
existence of some modulation mechanism in the space between the sun and the earth 
prior to impinging of the solar cosmic ray particles. In the November event, it seems 
that solar magnetic cloud by which the Forbush decreases were caused played an im- 
portant part as such a modulation mechanism. 


4, Concluding Remarks 


In the four unusual increase events of cosmic rays occurred in 1960, the increases 
observed in the polar regions showed the north-south asymmetry between both hemis- 
pheres, which property differs from event to event. It is suggested from these facts 
that polar distribution of the phenomenological aspects of the increases is determined 
by not only the isotropic solar cosmic ray particles but also the anisotropic particles 
deflected from their initial emitted direction from the sun. Some modulation on solar 
cosmic ray particles seems possible by considering the operation of the solar magnetic 
cloud by which the Forbush decrease was caused. 
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Appendix 


Data of cosmic ray neutron intensity obtained at Syowa Base on November 12 and 15, 1960 
are given in the following tables, where 15-minute values, corrected for pressure, scaling factor 8, 


are given. 


November 12, 1960 


UT 00—15 15—30 30—45 45—00 
11 94 93 97 92 
12 97 98 93 91 
13 96 93 103 98 
14 105 115 125 126 
15 137 134 161 159 
16 168 159 160 154 
17 158 159 158 161 
18 161 147 155 156 
19 137 bn = oy 
20 202 - 192 7 
21 180 172 160 164 
22 155 138 160 


November 15, 1960 


0 84 86 88 90 
1 82 83 89 82 
2 79 86 110 180 
3 160 159 164 157 
4 149 148 157 152 
5 146 142 148 140 
6 138 136 130 126 
7 130 118 118 118 


LETTER TO THE EDITORS 


Secondary Magnetization and the Palaeomagnetism 
of Some Pliocene Rocks of Japan 


(Received May 30, 1961) 


Abstract 


The effect of secondary magnetization on the apparent position and polarity 
of a palaeomagnetic pole is examined and applied to measurements made by 
Momose on Pliocene rocks. It appears that in the Lower-Middle Pliocene the 
pole was near 70°N, 150°W, and that the polarity in the Lower Pliocene was 
normal and in the Middle Pliocene was reversed. A reversal of the earth’s field 


appears to have occurred near the Lower-Middle Pliocene boundary. 
Introduction 


In studies of palaeomagnetism the positions of the earth’s magnetic poles in the 
geological past are computed from the observed directions of remanent magnetization 
in rocks, assuming a geocentric dipole field. The polarity of the field is called normal 
or reversed according to whether the inferred south-seeking pole lies in the northern 
or the southern hemisphere. Both polarities are commonly found in rocks of approxi- 
mately the same age from the same locality. 

It is necessary to establish that the observed direction of magnetization is the 
actual direction of the ancient field. One problem is to recognize and, if possible, to 
remove any secondary magnetization that the rock may have acquired in more recent 
times and which is superimposed on the original primary magnetization. This paper 
is concerned with the effects of secondary magnetization on the inferred positions and 
polarities of the ancient pole. We use the term secondary magnetization to denote 
a more or less stable component of magnetization whose direction and polarity are 
those of the earth’s average present field. The axis of this field coincides approxi- 
mately with the axis of rotation of the earth. It is not necessary here to specify the 
particular mechanism by which a secondary magnetization may be acquired but it 
should be noted that secondary magnetization may be magnetically hard and difficult 
to remove by demagnetization without simultaneously destroying the primary magnet- 
ization. It might perhaps be a chemical remanent magnetization or a viscous isother- 


mal magnetization resulting from the continual exposure of the rock to the earth’s 
field for many thousands of years. 


The Effect of Secondary Magnetization on Pole Position 


It is well known that secondary magnetization produces an apparent shift of the 
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palaeomagnetic pole along the meridian great circle that passes through the true 
palaeomagnetic pole. When the primary magnetization has reversed polarity and the 
secondary magnetization is large, the secondary magnetization may produce an apparent 
shift of the pole large enough to produce an apparent change of 180 degrees in the 
longitude of the pole, accompanied by an apparent change in polarity. 

The effects are most simply seen by considering the magnitudes and directions of 
dipoles at the centre of the earth which may be regarded as the sources of the fields 
which give rise to the primary, secondary and resultant magnetizations respectively 
in the rock at the place where it is situated. These dipoles are illustrated in Fig. 1. 


(a) (b) (c) 


Fig. 1. The effect of secondary magnetization on apparent pole position. The diagrams are 

drawn in the meridian plane containing the true palaeomagnetic pole P; 0-0 is an equatorial 

diameter. A represents the apparent pole in each case. In (a) the primary magnetization is 

formal; in (b) the primary magnetization is reversed and the secondary magnetization is small; 
in (c) the primary magnetization is reversed and the secondary magnetization is large. 


The diagrams of Fig. 1 are drawn in the meridian plane containing the dipole sources. 
_N and S are the present north and south geographic poles, C is the centre of the 
earth, and O-O is an equatorial diameter. The magnetic moment of the dipole corres- 
ponding to the primary magnetization is represented by the vector CX. The dipole 
moment corresponding to the secondary magnetization is represented by CY and is, 
of course, parallel to and has the same sense as the earth’s present magnetic axis. 
The resultant moment of CX and CY is represented by CZ. The points of intersection 
of the axes of CX, CY and CZ with the surface of the earth in the northern hemis- 
phere are the respective positions of the true palaeomagnetic pole P, the present geo- 
graphic north pole N, and the observed (or apparent) palaeomagnetic pole A. 

We consider points on the meridian NPOS that passes through the true palaeo- 
magnetic pole. An apparent pole observed to lie between N and P can only arise 
from the situation illustrated in Fig. 1(a), and must have normal polarity. An apparent 
pole between P and O, as in Fig. 1(b), must have reversed polarity. An apparent 
pole whose antipole in the southern hemisphere lies on the meridian NPOS, as in 


Fig. 1(c), must have normal polarity. 
Further, an apparent pole corresponding to Fig. l(a) can only arise from a true 


224 F.H. HIBBERD 


pole at P which had normal polarity. Apparent poles corresponding to Fig. 1(b) and 
1(c) can only arise from a true pole at P which had reversed polarity. 


Application to Some Pliocene Data 


We now use the above principles to examine some determinations by Momose (1958) 
from Pliocene volcanic rocks in Central Japan. The inferred pole positions are set out 
in stratigraphic sequence in Table 1, with normal polarity indicated by— and reversed 


Table 1. Pliocene pole positions, after Momose (1958) 
Komoro Group 


; Pole : 
Age | Symbol Locality one Long. Polarity 
1 Fukazawa 1 74 169W _ 
Lower Pliocene 2 Fukazawa 2 75 175E a 
a Oya 12 20E — 
4 Higashi zawa 12 44E : - 
5 Miyanoshita 8 171W “f 
Middle Pliocene 6 Sehaya gawa A 32 171W SC + 
7 Sehaya gawa 65 122wW | + 
8 Hanare yama 63 110W ao 
9 Yasuhara 67 121W + 
Shigarami Group 

Miocene a Konabe 26 119W _ 
Pape Oincant b Akefuji bashi 75 141W = 
c West of Okkayo 72 175W - 

d West of Akefuji 
bashi A 15 172E + 

e West of Akefuji 
ey 8 142W + 
Middle Pliocene f Kawashimo A 32 62E — 
g Kawashimo 17 50E a 
h Machi 15 157W + 
i East of Doai 35 142W oo 
Upper Pliocene j Kuroiwa 70 6W “+ 


polarity by+. They come from two groups of rocks, Komoro and Shigarami, which 
are believed to include the whole extent of the Pliocene. The rocks in the Komoro 
group show no sign of geological disturbance. Some samples in the Shigarami group 
were obtained from rocks that had been disturbed and corrections had to be made 
before computing the pole positions, but nevertheless the Shigarami results provide a 
useful check on the Komoro results. Each pole position is the average of several 
samples from each rock mass and the directions of magnetization within each set of 
samples in the Komoro group were concentrated within a range of 5°-7°. Within each 
group there is a considerable spread of the latitudes of the poles. Momose has shown 
that the longitudes of the poles lie roughly along a meridian great circle and he 
suggests that the positions of the poles may correspond to intermediate stages in a 
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reversal of the earth’s field. 

As an alternative interpretation, we shall show that the positions and polarities of 
the poles are consistent with a rapid reversal of the earth’s field which has left practi- 
cally no sign of intermediate stages, together with the superposition of a hard secondary 
magnetization in some specimens. Momose had subjected his specimens to A.C. 
demagnetization in a field of 400 oersteds so that any soft secondary magnetization 
must have been removed. However, there is still the possibility of a hard secondary 
magnetization being present. 

The pole positions in Table 1 are displayed by drawing the distributions in longi- 
tude and in latitude. Fig. 2 shows the longitudes of the pole positions in the north- 
ern hemisphere, in the nearest 10 degree intervals, and confirms that the pole posi- 
tions lie near a great circle passing through the present north pole. Most of the 
poles from both groups lie near the meridian 150°W. Poles 3, 4 and g, f, j lie near 
30°E in the northern hemisphere, so that their antipoles in the southern hemisphere 
lie néar 150°W. 

The latitudes of the poles that lie near the meridian 150°W, and the latitudes of 
the antipoles (in the southern hemisphere) of the poles that lie near the meridian 30°E, 
are shown in the nearest 10 degree intervals, in Fig. 3. The polarities of the poles 
are also shown. 

The polarities shown below the various latitudes at the bottom of Fig. 3 are 


Komoro 


Shigarami 


120 E 180 120 Ww 60 


Fig. 2. Longitudes of the northern hemisphere poles 

of the Komoro and Shigarami groups, showing that the 

poles lie near a meridian great circle passing through 
150°W and 30°E. 
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Komoro 


90 S 60 


Shigorami 


Fig. 3. Latitudes of the poles of Fig. 2 that lie near 
150°W and of the antipoles of the poles of Fig. 2 that lie 
near 30°E. The true palaeomagnetic pole appears 
to lie near 70°N, 150°W. 

those that the inferred poles should have in these respective latitudes if their positions 
are the result of a secondary magnetization superimposed on a primary magnetization 
corresponding to a pole of either polarity situated at 70°N. It is seen from Fig. 3 
and Fig. 1 that the positions and polarities of the poles of the Komoro group are 
completely consistent with the effects of secondary magnetization if the true magnetic 
pole was situated near 70°N, 150°W and if the polarity of the primary magnetization 
of specimens 1 and 2 was normal and that of specimens 3 to 9 was reversed. Except 
for poles a and j, which may have been affected by geological disturbance, the poles 
of the Shigarami group are also consistent with this interpretation if the primary 
magnetization of specimens b and c was normal and that of d to i was reversed. (In 
any case, pole j is younger than any of the others in either group). 

The fact that the inferred reversals of polarity between 2 and 3 of Komoro and 
between c and d of Shigarami appear to be simultaneous in the two stratigraphic 
sequences suggests very strongly that we have here an instance of a reversal of the 
field of the earth rather than a self-reversal of the directions of magnetization of the 
rocks. . 

We conclude that, if the presence of secondary magnetization is assumed, Momose’s 
observations are consistent with the earth’s magnetic pole in Lower-Middle Pliocene 
times being situated near 70°N, 150°W, that the earth’s field had normal polarity in 
Lower Pliocene times and reversed polarity in Middle Pliocene times, and that a rapid 
reversal of the field occurred near the Lower Pliocene—Middle Pliocene boundary. 


Reference 
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“ Geomagnetic Secular Variation during the Period from 1955 to 1960 2 TaNagata 
and Y. Syono (J. Geomag. Geoelect. 12, 84-98 (1961)). 
On 10th line from the bottom, p. 95, 
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